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Abstract 



The present paper deals with the classical trajectories for two super-integrable sys- 
tems : a system known in quantum chemistry as the Hartmann system and a system of 
potential use in quantum chemistry and nuclear physics. Both systems correspond to ring- 
shaped potentials. They admit two maximally super-integrable systems as limiting cases, 
viz, the isotropic harmonic oscillator system and the Coulomb-Kepler system in three di- 
mensions. The planarity of the trajectories is studied in a systematic way. In general, 
the trajectories are quasi-periodic rather than periodic. A constraint condition allows to 
pass from quasi-periodic motions to periodic ones. When written in a quantum mechanical 
context, this constraint condition leads to new accidental degeneracies for the two systems 
studied. 
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1. Introduction 



Three-dimensional potentials that are singular along curves have received a great deal 
of attention in the recent years. In particular, the coulombic ring-shaped potential Vq 
(see Eq. (3.1)), revived in quantum chemistry by Hartmann and collaborators [1], and the 
oscillatory ring-shaped potential Uq (see Eq. 2.1)), systematically studied by Quesne [2], 
have been investigated from a quantum mechanical viewpoint by using various approaches 
(including Schrodinger, Feynman path integral, nonbijective canonical transformations, 
and Lie-like approaches) . In this connection, we mention the works of Refs. [3-9] on the 
potential Vq. Furthermore, the potential Uq and various extensions of Uq and Vq have 
been worked out by means of algebraic, path integral, and Lie-like techniques [10-14]. 
Finally, the so-called ABC and ABO systems (which are companions of the Vq and Uq 
systems concerning an Aharonov-Bohm plus a Coulomb potential and an Aharonov-Bohm 
plus an oscillator potential, respectively) have been the object of several investigations in 
a nonrelativistic formulation [15-18] as well as in a relativistic one [19]. 

The Uq and Vq systems have been little examined from a classical viewpoint [7,20]. It 
is the aim of this paper to study the classical motion for a (charged) particle embedded in 
the potential Uq or Vq. More specifically, we shall investigate the planarity, the periodicity 
and the semi-classical quantization of the bounded motions for the Uq and Vq systems. 

Among the potentials in n = 3 dimensions, the potentials Uq and Vq have a spe- 
cial status in the sense that they correspond to super-integrable hamiltonian systems with 
n + 1 = 4 integrals of motion. The latter two systems occupy a position intermediate 
between integrable hamiltonian systems (with n = 3 constants) and maximally super- 
integrable hamiltonian systems (with 2n — 1 = 5 constants). In addition, when the param- 
eter Q vanishes, the Uq and Vq systems become maximally super-integrable systems ; as a 
matter of fact, Uq = q corresponds to the oscillator system and Vq=q to the Coulomb-Kepler 
system. Let us recall that, in n dimensions (n arbitrary), four maximally super-integrable 
hamiltonian systems are known : the oscillator system, the Coulomb-Kepler system, the 
Calogero-Moser system, and a system introduced by Smorodinsky and co-workers [21,22] 
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(see also Ref. [23]). 



The Uq and Vq systems thus turn out to be two interesting laboratories among the 
super-integrable systems. As ring-shaped systems, they may play an important role in all 
situations where axial symmetry is relevant. For example, the Vq system is of interest for 
ring-shaped molecules like cyclic polyenes [1]. Further, the Uq system is of potential use 
in the study of (super-) deformed nuclei. 

The plan of the present paper is as follows. For each of the potentials Uq and Vq we 
study in Sections 2 and 3, respectively, the constants of motion, the classical trajectories, 
their planarity, (quasi-) periodicity and quantization. The quasi-periodicity of the bounded 
motions is connected in Section 4 to the concept of "local" symmetry of a spectral problem 
as first introduced by Coulson's school in connection with the Kepler problem [24]. The 
planarity of the bounded motions is studied with the help of some new formulas, rele- 
gated to the appendix, for the torsion (and the curvature) of the trajectories of a given 
hamiltonian system. 



2. The Oscillatory Ring-Shaped Potential 



2.1. Generalities 



We shall deal in this section with the potential (energy) 

U Q = ln 2 (x± 2 + x 2 2 + x 3 2 ) + Iq^- j O>0 Q>0 (2.1) 

2 2 xi z + x 2 

which clearly exhibits an 0(2) cylindrical symmetry. The limiting case Q = corresponds 
to an isotropic harmonic oscillator (in R 3 ) and will serve for testing the results to be 
obtained. 

The potential Uq is a special case of the potential 

V 3 = a ( Xl 2 + x 2 2 + x 3 2 ) +( 3 -L + h(^) 2 1 2 (2.2) 

X 3 Z Xi Xi 2j rX 2 2 



of Makarov et al. [22] . The potential V 3 allows the separation of variables for the Hamilton- 
Jacobi equation and, also, the Schrodinger equation in four systems of coordinates, viz, 
the spherical, circular cylindrical, prolate spheroidal and oblate spheroidal coordinates. 
According to Makarov et al. [22] , there are four functionally independent integrals of motion 
for V 3 and, thus, for Uq. More precisely, the four integrals for Uq can be constructed in 
terms of the relevant Hamilton function and the A^s listed below [2,20,22]. 

(i) Spherical coordinates : 

x\ = r sin#cos</? x 2 = r sin6>sin</? x 3 = rcos9 

(2.3a) 

a 1 = e, 2 + e 2 2 + e 3 2 + Q — L- a 3 = e 3 . 

sin 



(ii) Circular cylindrical coordinates : 

xi=pcosnp X2 = psimp x 3 = z 

(2.3b) 

A 2 = \ ( P 3 2 + ftW) A 3 =£ 3 . 



(iii) Prolate spheroidal coordinates : 

x\ = a sinhry sin a cos <p x 2 = a sinh ?y sin a sin x 3 = a cosh rj cos a 

(2.3c) 

A 4 = i x 2 + £ 2 2 + i 3 2 - a 2 ( Pl 2 + p 2 2 ) - a 4 2 sinh 2 V sin 2 a + Q ^ I 7 ~ - A 3 = i 3 . 

sinh rysin a 



(iv) Oblate spheroidal coordinates : 

xi = a cosh ?y sin a cos </? x 2 = a cosh 77 sin a sin <p x 3 = a sinhry cos a 



5 



(2.3d) 

A, = £ 2 + £ 2 + £ 2 + a 2 ( Pl 2 + p 2 2 ) + a 4 Q 2 cosh 2 V sin 2 a + Q cosh " ^ + si f a a 3 = ^ 

cosh 77 sin a 

(It is important to note that, in order to adhere to standard notations, terms of type (j>k) 
and (— Mfc or Lk) in the notation of Makarov et al. [22] or Kibler and Winternitz [7] are 
replaced here by (ipk) and (— i£k), respectively. The variables pk and £k in the present 
paper stand for linear and angular momenta, respectively, in atomic units.) 

Returning to Cartesian coordinates, we can write the first three integrals of motion 

as 

A (i 2 , n 2 , n 2 , r\ X ^ + + Xz A 1 / 2 , r>2 2\ a 

A 1 =£ 1 +£ 2 +£ 3 +Q 5- 5 — A 2 = - {P3 +11x3) A 3 = xip 2 -x 2 pi. 

xi 2 + x 2 2 

(2.4) 

The remaining two then satisfy 

A 4 = A 1 + 2a 2 (A 2 - H) A 5 = A 1 - 2a 2 (A 2 - H) (2.5) 

with 

H=\{p 1 2 +p 2 2 +P3 2 )+U Q . (2.6) 

We see that {H, Ai, A 2 , A3} form an integrity basis for all integrals of the considered 
Hamiltonian system. 

A consequence of Eq. (2.5) is that any potential that allows the separation of vari- 
ables in spherical and cylindrical coordinates will also allow the separation in prolate and 
oblate spheroidal coordinates. This holds both for the Hamilton-Jacobi and Schrodinger 
equations. 

2.2. Equipotential Surfaces 

The equipotentials corresponding to given values of Uq, with Uq > O y^, are de- 
scribed (in circular cylindrical coordinates) by 

Z = ± h,~p V" 02 P 4 + 2U QP 2 -Q p=Vxi 2 +%2 2 z = x 3 (2.7) 
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with the restrictions 

P< <P 2 < P> (2.8) 

where 

2 _ u Q - \Ju q 2 -wq 2 _ u Q + ^Ju q 2 -wq 

P< ^ P> ^ • (2-9) 

All the equipotentials corresponding to bounded values of Uq with Uq > y/Q are 
bounded surfaces which collapse into simple lines (actually circles in the z = plane) 
for Uq = O y/Q. Of course, each equipotential reduces to a simple sphere, of radius 
y/2 Uq/Q, in the limiting situation (Q = 0, Uq > 0). 



2. 3. Trajectories 



We now derive the equations of motion for a particle with potential energy Uq. We 
use here the cylindrical coordinates q = (p, <p, z) in which the Hamilton function H for a 
particle of (reduced) mass \x = 1 in the potential Uq reads 

H = H(p, q)= 1 - (p 2 + 1 p 2 + p 2 ^ + \ fi 2 (p 2 + z 2 ) + \ Q 1 (2.10) 

where p = (p p = p, p<p = P 2( Pi Pz = are the relevant canonically conjugated momenta. 
The corresponding Hamilton- Jacobi equation can be written 



dS_ 1 
~dt + 2 



dp 



+ 



1 /<9S 



+ ^ 2 ( p 2 + , 2 ) + 1q-1 = o. 



(2.11) 



Obviously, <p is a cyclic variable so that we put 

dS 



p v = tt- = m to G R 

dtp 



(2.12) 



and, since the system of the particle in Uq is conservative, we look for a solution of (2.11) 
in the form 

S = Si(p) + mtp + S 2 (z) - Et. (2.13) 
This leads to the separated equations 



1 fdSiY 1 2 2 1. 2 ^1 



Kt) ^^ 2 =— (-) 



which are coupled via a (positive) separation constant K. Equations (2.14) show that the 
only possible motions are bounded motions. They occur for 



E > K > K 2 - tt 2 M 2 > 



(2.15) 



where 



Indeed, Eqs. (2.15) ensure that 



M 2 = m 2 + Q. 



-z < z(t) < z p x < p(t) < p 2 



(2.16) 



(2.17) 



with 



>J1{E ~ K) 



Pi 



92 



Vk + VW- fi 2 M 2 



(2.18) 



Direct integration of (2.14) yields 

Slip) = \n vV - pi 2 ) (P2 2 - p 2 ) 



1 K 

+ 2 n sin " 



2p 2 - (Pi 2 + P2 2 ) 

P2 2 ~ Pi 2 



- - \M\ sin" 1 
2 1 1 



(pi 2 + P2 2 ) p 2 -2 Pl 2 p 2 < 



(P2 2 ~ P1 2 )P 2 



S2(Z) 



E-K 



: 'i-m +sm- i/ : 



(2.19a) 
(2.196) 



Equation (2.19a) is valid for K 2 - tt 2 M 2 > ; in the special case K 2 - Vl 2 M 2 = 0, we 
have 

^[k \m\ Um 

=>- Si(p) = constant 



Pl = „ 2 = „„ = — = 7 = 

and the motion is constrained to p = pq. 



(2.20) 



The action S follows from (2.13) and (2.19). Therefore, the equations of motion 

(2.21) 



9S -8 



— =(32 —=(3s 
dK dm 



can be derived in a straightforward way. In fact, in view of the relations 



. dS 2 dS, dS 2 dSi dS 2 

Pi = T77T ~t P2 = + 



dE 



dK dK dK dE 



dm 



(2.22) 



we choose t' Q = —fii, t = —(Pi + /3 2 ) and <p = 03 for constants of motion. Thus, we 
obtain (for K 2 - 2 M 2 > 0) 



i_to= 2S] Sin " 



1 m 



2 M 



sin 



2p 2 -(pi 2 + p 2 2 ) 

P2 2 - Pi 2 



(pi 2 + P2 2 ) p 2 -2pi 2 pS 

(P2 2 ~P1 2 )P 2 



(2.23) 



from which we can extract the coordinates p(t), (fit) and z(t). We finally arrive at 
p(t) = ±= a/Pi 2 + P2 2 + (P2 2 - Pi 2 ) sin[20(t - t )] 



<p(t) = ipQ + 



1 m . _x f p 2 2 ~ Pi 2 + (pi 2 + P2 2 ) sin[2Q(t - t Q )] \ 

2 |M| Sm 1 p! 2 + p 2 2 + (p 2 2 - p! 2 ) sin[20(t - t )] J 

= 2 sin[O(t-t )] 



(2.24a) 

(2.246) 
(2.24c) 



in agreement with a previous derivation based on functionally independent constants of 
motion [20]. Equations (2.24a) and (2.246) are valid for K 2 - tt 2 M 2 > and also for 
K 2 — Q 2 M 2 = under the condition pi = p 2 = po- Note that in (2.246), the parameter 
m/\M\ (which is ±1 in the limiting case Q = 0) is given by 

,2 



1 



M 2 Q 2 Pl 2 p 2 2 

so that the variables p, <p and z depend on the six constants p±, p 2 , <po> z o, to and t' Q . 



(2.25) 



At this stage, it is worth mentioning that (2.246) is equivalent to 

1 m m _! ( (p! 2 + p 2 2 ) tan[Q(t - t )] + p 2 2 - pi' 

(pit) = cp 1 r a + 1 r tan < 

rK 1 ru 2 |M| |M| \ 2pip 2 



with 



sin a 



2 2 
P2 ~Pi 

Pi 2 + P2 2 



cos a 



2piP2 



(2.26) 



(2.27) 



Pi 2 + P2 2 ' 

Indeed, Eq. (2.26) naturally arises when solving the Newton equations for the dynamical 
system under consideration. 
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As a resume, all the trajectories are bounded. The most general motion corresponds 

to 

O>0 Q>0 E>K>n\M\>0. (2.28) 

It follows from (2.17) that the classical trajectories lay between two cylinders of height 2z 
and radii p\ and p 2 . The real constant m is the ^-component £3 of the angular momentum 
for the system being studied ; we thus recover that M 2 = A 3 2 + Q is a constant of motion 
(cf. Ref. [22]). Furthemore, the positive constant K is another constant of motion, namely, 
K = E — A 2 . In terms of pi, p 2 and zo, the three constants of motion m, K and E are 
given by 

m 2 = O 2 Pl 2 p 2 2 -Q K=^Q 2 ( Pl 2 + p 2 2 ) E=^Q 2 ( Pl 2 + p 2 2 + z 2 ). (2.29) 

Finally, it is to be emphasized that the trajectories are neither periodic nor planar in 
general. 

2-4- Planarity of Trajectories 

The question of planarity of (bounded and nonbounded) motions can be tackled by 
looking at the torsion r of the trajectories. A planar motion corresponds to r = along 
the trajectory. For the Uq system, by using the formula for r = — NUM/DEN derived in 
the appendix, we get 

p 8 x NUM = Qm[(Q- Q 2 p 4 )z + 4 tt 2 p 2 (xx + yy)z] 

p 8 x DEN = (Q -fi 2 p 4 ) 2 (m 2 + p 2 z 2 ) (2.30) 
+ 2 fiV (Q - OV )(££ + yy)zz + (0 2 p 4 ) 2 (x 2 + y 2 )z 2 . 

The condition r = can be seen to be equivalent to 



Qm[(Q- Q 2 p 4 )z + 4 Q 2 p 3 pz] = Q 
with 



in 



(Q-n 2 P *)f t (z)-zf t (Q-n 2 p*) 



(2.31) 



pp = xx + yy = ±y / [2E - z 2 - fl 2 (p 2 + z 2 )] p 2 - M 2 . (2.32) 
Solutions of Eq. (2.31) are : 
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(i) (5 = 0. In this case the studied system reduces to the harmonic oscillator for which 
of course all trajectories are planar (ellipses with the attractive "sun" in the center) . 



(ii) z = z = 0. We then have z = in (2.24c) and hence E = K, see (2.29). 

(iii) m = 0. The motion is restricted to the plane ip = </?o, see (2.24b) and (2.26), and 
the limiting radii satisfy p\p 2 = y/Q/Q,, see (2.29). 

(iv) p 4 = Q/Q 2 . Indeed, the generic solution of (2.31) is Q — 2 p 4 = cqz with 
Co = constant. However, since z and p have different periods, see (2.24), the latter solution 
can only hold for p = pi = P2 = (Q/Q 2 ) 1 ^ 4 which implies m = 0, see (2.29), and hence 
either zq = (no motion) or cq = 0. 

To sum up, excluding the well known harmonic oscillator corresponding to the limiting 
case Q = 0, we find that the trajectories are planar only if we have m = or E = K. 

2.5. Periodicity and Quasi- Periodicity 

From Eq. (2.23), we easily see that the projections onto the xy plane of the trajectories 
are described by 

p = V2 P1P2 (2.33) 

'pi 2 + P2 2 - {p2 2 - Pi 2 ) sin 2 ^ (<p - <fi ) 



and, thus, these projections are not closed in general. Equation (2.33) shows that the 
projection onto the xy plane of a given trajectory is closed if the condition 

\Ml = h kieZ k 2 eZ (2.34) 
m K2 

(with k\ 2 > k 2 2 ) is fulfilled. The situation ki 2 = k 2 2 = I corresponds to the limiting case 
Q = for which we know that all trajectories (actually ellipses) are planar and periodic. 
In the situation where (2.34) is satisfied, then not only the projection, onto the xy plane, 
of the trajectory is periodic, but the trajectory itself is periodic too. The period of the 
motion is then 

T = k 1 T T o = ~^ (2-35) 
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where To is the oscillator period corresponding to the limiting case Q = 0. (See also 
Ref. [20] for an alternative derivation of (2.35).) 



In the general case Q > 0, when the condition (2.34) is satisfied, we have 

2 r\ i 2 



m 2 = Q 1 2 " z 2 piW = ^ g" 1 2 (2.36) 



_&2 2 2 _ Q_ fa' 

h 2 - k 2 2 Pl 92 ~ O 2 kl 2 - k 2 

for the periodic motions. The requirement that the trajectories be periodic thus leads 
to a "quantization condition" for the component £3 of the angular momentum and the 
mean geometric radius v /plp2. Motions for which (2.34) is not satisfied are quasi-periodic 
motions rather than periodic ones. 



The potential energy 



Q 



1 r To 

Jo 



dt (2.37) 



averaged over a duration To (a common period for p and z) can be easily calculated from 
(2.24a) and (2.24c). We get 

U Q ]=\si 2 {p 1 2 + p 2 2 + zo 2 ) + \Q— => Uq:= 1 -E + ]-Q— (2.38) 

4 I pip 2 I I P1P2 

for periodic and quasi-periodic motions. Thus, the virial theorem (T = V = E/2) for the 
three-dimensional isotropic harmonic oscillator applies to the Uq system only when Q = 0. 

2.6. Particular Cases 

The case (5 = 0. As a general check of the correctness of (2.24a), (2.246) (or (2.26)) 
and (2.24c), we can verify, by calculating x = pcostp and y = p sin 99, that the limiting 
case Q = effectively corresponds to a three-dimensional isotropic harmonic oscillator with 
angular frequency O. As a further check, it can be verified that for Q = the components 
of the angular momentum for the considered particle are 

£2 = - z a/pi 2 + P2 2 |cos(v? + 77) cos[0(t[, -t )] - sm(tp - ^) sin[0(to - £ )]} 
£3 = m 

li = +-^=z a/pi 2 + p2 2 {sin(</? + ^) cos[0(to - t )] + cos(y? - ^) sin[0(to - t )] | 

(2.39) 
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indicating that the elliptic trajectory corresponding to pi, p 2 , <fo, zq, to and t' is in the 
plane perpendicular to the constant vector (^1,^2,^3)- 



The case K = E. Equations (2.18) and (2.24c) show that the trajectories are in the 
xy plane when E = K, a situation that corresponds to case (ii) of section 2.4. They are 
described by (2.33) and thus are not closed in general. However, the planar trajectories 
corresponding to E — K become closed if the condition (2.34) is satisfied. 

The case m = 0. Case (iii) of section 2.4 indicates that the trajectories are planar 
for m = and, from Eqs. (2.24), we see that they are in the pz plane corresponding to 
(p(t) = (fQ. All the trajectories (p(t), z(t)) are periodic, of period To- 

The case K = 0|M|. The two cylinders of radii p\ and P2 collapse into a single one, 
of radius po, when K = Q\M\ ; therefore, the trajectories (neither closed nor planar in 
general) are on a cylinder of height 2z and radius p (see (2.17), (2.18) and (2.20)). Then, 
Eqs. (2.24) lead to the Cartesian coordinates 

( 777- \ / Tfl \ 

(2.40) 

where ip is some new constant. The motion thus results from the combination of a circular 
motion of angular frequency 

^ n W\ = ^ <2M) 

in the xy plane and of an oscillatory motion of angular frequency O in the z direction. 
(Without loss of generality, we can take m > ; the case m = is trivial.) In the limiting 
case Q = 0, the trajectory reduces to an ellipse, the projection of which on the xy plane is 
a circle of radius pq. 

In the case K = Q\M\ and Q > 0, the trajectories are not closed in general. However, 
if we further assume that the condition (2.34) is satisfied, the trajectories become closed 
and of period T = k\To- They are given by (2.40) and (2.41) with 



n _ h 

u) k 2 




(2.42) 
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with fei > k 2 , h G N - {0, 1} and fc 2 G N — {0}. 



2.7. Semi- Classical Quantization 

We now use the Bohr-Sommerfeld-Kramers quantization conditions to derive, in a 
semi-classical way, the quantum mechanical spectrum for the Uq system. Here, these 
conditions read 



j> Ppdp = (n p + e p )2n n p eN j) p z dz = (n z + e z )2n n z eN 



(2.43) 



where e p and e 2 stand for (rational) numbers to reconcile the old and new theories of 
quanta. (In the case of a one-dimensional harmonic oscillator, we know that e z = 1/2.) In 
other words, we have 

2 [SM - Si(pi)] = (n„ + e p )2n 2 [S 2 (z ) - S 2 (-z )] = (n z + e z )2n (2.44) 

and, from (2.19), we obtain 



1K !|^: E-K 



m G Z n p GN n 2 GN 



Consequently, the eigenvalue E is given by 

E= (\M\+2n p + n z + Cl \M\ = y/m 2 + Q 

(2.46) 

in agreement with the result by Quesne [2]. In fact, Eq. (2.46) demands that e z +2e p = 3/2, 
a relation which is reminiscent of the three-dimensional isotropic harmonic oscillator. 

3. The Coulombic Ring-Shaped Potential 



3.1. Generalities 



The coulombic ring-shaped, or Hartmann, potential (energy) is 

V Q = -Z^= 1 +\Q 2 ] Z>0 Q>0 (3.1) 

Vx{ 2 +X 2 2 +X 3 2 2 X! 2 + X 2 2 

14 



where Z = i]a 2 and Q = q rf a 2 in the notation of Hartmann [1] and of Kibler and Negadi 
[3]. Such an 0(2) invariant potential reduces to an attractive Coulomb potential (in R 3 ) 
in the limiting case Q = and this will prove useful for checking purposes. 

Clearly, Vq is a special case of the potential (in spherical coordinates) 

1 cos# 1 . . 

r r l sin r 1 sin 

introduced by Makarov et al. [22]. Therefore, the Schrodinger equation and, thus, the 
Hamilton-Jacobi equation for Vq are separable in spherical and parabolic rotational coor- 
dinates. In this respect, let us simply recall that there are four functionally independent 
integrals of motion for Vq which can be obtained from the relevant Hamilton function and 
the Sfc's given below [7,20,22]. 

(i) Spherical coordinates : 

x\ = r sin#cos</? X2 = r sin(9sin</? X3 = rcos6* 



(3.3) 



b x = e x 2 + e 2 2 + h 2 + Q -Aa- B 2 = h 

sin^ 6 



(ii) Parabolic rotational coordinates 



x\ = Vab cos ip X2 = Vab sin ip x^ = — (a — b) 



b 3 = 1 (hP2 + P2h - hPi - Pih + 2 z ^-J - Q (L ~r ) B ' = ' 

2 V a + b ab 



(3.4) 



(The constant B% is given in a quantum mechanical form. The factor +2[. ..] in the 
corresponding expression of Ref. [22] for B s should be changed into — 2[. . .].) 



3.2. Equipotential Surfaces 
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The equipotentials are formally given (in circular cylindrical coordinates) by 



z = ± 



y^ZV - (Q ~ 2V QP 2 ) 2 p= v / ^i 2 + ^2 2 z = x 3 . (3.5) 



Q - Wqp 1 

In contradistinction with Uq, three cases must be considered for the potential Vq. 

(i) The case — (1/2)(Z 2 /Q) < Vq < : The equipotentials are given by (3.5) with the 
restrictions 

. . g ~ V^ 2 + 2QV Q g + y^ 2 + 2QV Q 
P<<P<P> P< = P> = (3-6) 

and each equipotential is a bounded surface which reduces to a circle in the z = plane for 
Vq = —(1/2)(Z 2 /Q). In the limiting case Q = 0, each equipotential for V"o < becomes a 
sphere of radius —Z/Vq. 

(ii) The case Vq = : The equipotentials are obtained from 

^ = ±lp v / 4^V-Q 2 (3.7) 



with the restriction 



P >\% (3.8) 



and thus the equipotentials are not bounded surfaces in this case. (The limiting case Q = 
yields of course a sphere of infinite radius.) 

(iii) The case Vq > : The equipotentials are given by (3.5) with the restriction 



, „ -Z + + 2QVq 1 Q .... 

p<<p<p> p< = ^ p> = viV^ ( > 

and here again all the equipotentials are bounded surfaces. 
3. 3. Trajectories 

3.3.1. Preliminaries. We shall solve the Hamilton- Jacobi equation for Vq in spherical 
coordinates. The general pattern to be followed for Vq resembles the one for Uq so that 
we shall mention only the basic steps. 
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The Hamilton-Jacobi equation for a particle of reduced mass \i = 1 in Vq can be 
written as 



dS_ 1 
~dt + 2 



\dr J 



89 



+ 



1 fdS 



- 



* 2 sin 9 \d(p 



r 2 r 2 sm 2 6> 



(3.10) 



and a solution of the type 



S = S 1 (r) + S 2 (9) + mp-Et 



(3.11) 



is easily found by solving the two equations 

2 



dSi 
dr 



2Er 2 - 2Zr = -K 



- 1 + (m 2 + Q)— — = K 



de 



sin 2 6 



(3.12) 



where iiT is a (positive) separation constant. Equations (3.12) show that two kinds of 
motions may occur here : nonbounded motions for E > and bounded motions for 



1 Z 2 

2 ~K 



< E < 0. 



(3.13) 



We shall restrict ourselves to the bounded motions and to the separatrix which corresponds 
to E = 0. As in the case for Uq, the equations of motion are given by (2.21) and we shall 
use the constant M defined via (2.16). 



3.3.2. Bounded motions. They are confined in a region delimited by 



Z - V^ 2 + 2EK Z + VZ 2 + 2EK 

fl < r \t) <T2 T\ = — T2 



-2E 



-2E 



9o < 6(t) <tv-9 



sin 6*o 



IMI 



K 



By assuming that -(1/2) (Z 2 /^) < E < 0, Eqs. (3.12) admit the solutions 



S 1 (r) = (-2E) l l 2 ^(r- ri )(r 2 -r) 



+ (-2E)- 1 / 2 Z sin -1 



2r - (ri + r 2 ) 
r 2 - ri 



- v^sin -1 



(ri + r 2 )r - 2rir 2 
(r 2 - n)r 



(3.14) 



(3.15a) 
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S 2W JM i sin- 



1 



- VK sin -1 



cos 6*o 
cos 9 
cos 6>o 



-1 + 



sin 2 6>o 
1 — cos 9 



— sin 



cos 6>o 



-1 + 



sin 2 #o 
1 + cos 9 



Putting t = — Pi, Po = —2\fKf3 2 and (fo = (3s, we get from (2.21] 



t-t = -{-2E)- l/2 ^/{r - n)(r 2 - r) + Z(-2£)" 3/2 sin" 1 



2r - (n + r 2 ) 



r 2 - ri 



. _ i / cos 9 v . _ | 
Po = sin ( — ) + sin 



<p-<p 



1 m 



2 M 



sin 



cos 6*0 



cos 6>o 



-1 + 



(ri + r 2 )r - 2nr 2 
(r 2 - n)r 



(3.156) 

(3.16a) 
(3.166) 



sin 2 9q 
1 + cos 9 



— sin 



cos 9 



-1 + 



sin 2 6>o 
1 — cos 9 



(3. 



6c) 



Equation (3.166) is amenable to the form 

rcos9 = — |[2rir 2 — (r\ + r 2 )r] cos/3o + 2Jr\T2 \/ {r — r\) (r 2 — r) sin (3q } (3.17) 

r 2 — n L * J 

and, therefore, Eqs. (3.16) are in accordance with the ones derived by Kibler and Winternitz 

[20] from functionally independent constants of motion. Equations (3.16) are valid for 

-(1/2)(Z 2 /K) < E < 0. The case E = -(1/2)(Z 2 /K) < deserves a particular study 

(see section 3.6). In Eq. (3.16c), note that the parameter m/\M\ (which is ±1 in the 

limiting case Q = 0) can be obtained from 



m 



1 Q n + r 2 1 



M 2 " 2 Z nr 2 sin 2 9 
so that the variables r, 9 and <p depend on the six constants r±, r 2 , 6*o, (fo, to and (3q. 



(3.18) 



To sum up, the most general bounded motion corresponds to 

1 Z 2 



Z>0 



Q>0 



2 K 



< E < 



(3.19) 



and all finite trajectories take place between two spheres of radii T\ and r 2 (see (3.14)). 
The interpretation of the constant m (or M) is similar to the one for Uq : we have 
M 2 = B 2 2 + Q. Furthermore, it can be verified that K = B\. In terms of ri, r 2 and ^o, 
the three constants of motion m, K and E can be deduced from 



.2 nr 7 J^ sin 2 9o _ Q 

ri + r 2 



m 2 = 2 Z 



K = 2Z 



nr 2 
n + r 2 



E = —Z 



ri +r 2 



(3.20) 
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Here again, the finite trajectories are neither periodic nor planar in general. 



3.3.3. The separatrix. By taking E = 0, Eqs. (3.12) lead to 

1 K \M\ 
r(t)>r min r min =-- 9 < 6(t) < it - 9 sm9 = (3.21) 



so that the motions are nonbounded. From (3.11) and (3.12), it appears that the coordinate 
<p{t) is still given by (3.16c) when E = 0. In addition, the coordinate r(t) can be determined 
from the cubic equation 

r 3 + 3 r 2 r min - 4 r min 3 - (t - t ) 2 = (3.22) 

while the coordinate 0(t) follows from 

00 = sin" 1 ( + 2 tan" 1 ( J — - l) (3.23) 

\COs6oJ VV r min / 

or alternatively 

r cos 9 = cos 6 [(2 r min - r) sin (3 - 2 v/r~ Vr - ^min cos /3 ] (3.24) 
which parallels Eq. (3.17). 

3.4- Planarity of Trajectories 

The torsion r for the Vq system can be calculated from the general formula r = 
— NUM/DEN given in the appendix. We thus get 



p 8 x NUM = Qm{(Q- Z^j z + Z^ [(xx + yy) (r 2 + 3^ 2 ) - 3zzp' 



p 8 x DEN =\Q- Z'-^j (m 1 + p z z 2 ) (3.25) 

+ 2Z^ (q - Z^j (xx + yy)zz + (z^j + 2/ V- 
The condition r = for planarity amounts to solve 
Q m 

\ r" j r ' 

,P 4 \ d O] ( ryP^ 



Q-z£)z+z£(4pr-Srp)z 



Q m 
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(3.26) 



= 



with 



pr{Apr — 3rp) = (xx + yy)(r 2 + 3z 2 ) — 3zzp 2 



pp = xx + yy = ±J 2E - i 2 + 2Z- p 2 - M 2 . 




(3.27) 



(In (3.25) - (3.27), r and p refer to spherical and cylindrical coordinates, respectively.) 
Equation (3.26) is satisfied for all t in each of the following cases : 

(i) Q = 0. This is the case of the Kepler planetary system in which of course all 
trajectories are planar, namely, hyperbolas, or a parabola, or ellipses with the sun in one 
focus. 

(ii) z = i = 0. The integrals of motion satisfy \M\ = \[K and 6q = tt/2. 

(iii) m = 0. The motion is restricted to the plane cp = <po and we have \M\ = \fQ. 

(iv) p 4 = (Q/Z)r 3 . Indeed, case (iv) implies that m = (at least for bounded 
motions). Furthermore, the generic solution of (3.26) requires that Q — Z(p 4: /r 3 ) = c$z 
with Co = constant, a relation which does not hold for t arbitrary (except for the static 
case corresponding to z = and p 4 = (Q/Z)r s ). 

Consequently, in addition to the limiting case Q = (case (i)), the bounded motions 
are planar only if \M\ = \[K (case (ii)) or m = (case (iii)). 

3.5. Periodicity and Quasi- Periodicity 

In the case of finite motions, Eq. (3.16a) allows to obtain r as function of t by means 
of a transcendental equation. The latter equation indicates that r{i) is a periodic function 
of period 



a result in agreement with the fact that the parabolic coordinates a(t) and b(t) (with 
r = (a + 6)/2) are also of period Tc [7]. It should be noted that the period Tc bears 
exactly the same form as the one for the Coulomb-Kepler problem that corresponds to the 
limiting case Q = 0. Further, from (3.166) or (3.17) it is possible to show that cos[0(£)] (and 



T c = 2ttZ (-2E)~ S/2 



(3.28) 
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therefore sin[0(t)]) is a periodic function of period Tq. This result agrees with our previous 
work [7] due to the passage formulas cos6> = (a — b)/(a + b) and sin6> = 2\fab/{a + b). 
Finally, the consideration of (3.16c) leads to 

<p{t + hT c ) = <p{t) + 2 h 7T — h e Z (3.29) 

and, therefore, the (global) motion is periodic of period 

T = k l T c (3.30) 

if the "quantization" condition (2.34) for \M\/m is fulfilled (see also Ref. [20]). As a check, 
in the limiting case Q = 0, the global period T of the motion is simply Tq- 

In the general case Q > 0, when the condition (2.34) is satisfied, we have the following 
"quantized" expressions 



m =Q — 7) 7T sin Vo = - — — o o (3.31) 

h 2 -k 2 2 n+r 2 2Zk 1 2 -k 2 2 



for the periodic bounded motions. Here again, the bounded motions for which the condition 
(2.34) is not satisfied are quasi-periodic motions. 

The potential energy Vq averaged over the period Tc (of r and sin 9), namely, 



cannot be calculated easily from (3.16a) and (3.166). However, we can calculate in a 
straightforward way the virial of the Vq system. This leads to 

1 1 C Tc df 111 

V Q = 2E+-Q— -^- = 2E + Q^—^=—— (3.33) 

2 T c Jo r 2 sin 2 9 r x + r 2 ^/r x r 2 sin 9 

for periodic and quasi-periodic bounded motions. Therefore, the virial theorem (T = 
— V/2 = —E) for the three-dimensional Coulomb-Kepler problem applies to the Vq system 
only when Q = 0. 

3.6. Particular Cases 
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The case K = M 2 . Equations (3.12) show that the trajectories are in the xy plane 
when \M\ = \fK, a situation that corresponds to case (ii) of section 3.4. The finite 
trajectories are then described by 



p = 2 



Plp2 



Pi + P2 ~ (p2 ~ Pi) sin ^-((f- V?o)] 



(3.34) 



where p\ and p 2 are given by (3.14) with r = p. These finite trajectories are closed only if 
the condition (2.34) is satisfied. On the other hand, the separatrix corresponds to 



P = Pmin Si + tan 



where p m i n is given by (3.21) with r = p. 



1 \M\ 



2 m 



O - <Po) 



(3.35) 



The case m = 0. The trajectories are planar for m = (cf. case (iii) of section 3.4). 
They can be obtained in the plane ip{t) = tpo from (3.16a) and (3.17) and are all periodic, 
of period Tq. 



The case 2EK = —Z 2 . The finite trajectories are restricted to a sphere of radius 



K 



n = r 2 = r 



Z -2E 



K 
-2E 



(3.36) 



when E = —(1/2)(Z 2 /K). These trajectories (not closed in general) are described by 



r = r z = r cos 6> cos ( 2-k^- + V>o 

J-c 



m 

W\ 



tan 



tan 



(2tt^ + Vo 
sin 6>o 



(3.37) 



where ip denotes a new constant. They are periodic only if the condition (2.34) for \M\/m 
is satisfied. 

3.7. Semi- Classical Quantization 



The Bohr-Sommerfeld-Kramers quantization conditions 



p r dr = (n r + e r )2iv n r e N 



ped9 = (ne + eg)27i ng E N 



(3.38) 
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applied to the Vq system lead to 

2 [Si(r 2 ) - S 1 (r 1 )] = (n r + e r )2n 2 [S 2 (n - ) - S 2 (0 O )] = + e,)2yr. (3.39) 
By using Eqs. (3.15), we get 

(-2E)~ 1/2 Z- \[K = n r + e r - \M\ + Vk = n e + e e (3.40) 

from which we recover the eigenvalue 

1 pp. 

E = ---—. — |M| = ^m 2 + Q meZ n r G N n e EN 

2 (|M| +n r + n e + l) 2 v 

(3.41) 

as first derived by Hartmann [1]. Equation (3.41) implies that we must take eg + e r = 1, 
a relation which reflects the occurrence of a pair of two-dimensional harmonic oscillators 
in the treatment of the Coulomb-Kepler problem by means of the Kustaanheimo-Stiefel 
transformation [3,7]. 

4. Closing Remarks 



We have concentrated in this paper on the classical motions for two super-integrable 
systems, viz, the Uq and Vq systems. The motions are always bounded for Uq and can be 
bounded or nonbounded for Vq just as in the limiting case Q = that corresponds to the 
oscillator and the Coulomb-Kepler systems in three dimensions. The bounded trajectories 
are confined inside simple surfaces : cylinders for Uq and spheres for Vq. Both for Uq and 
Vq, the trajectories are planar only in a few situations : (i) naturally, when (5 = 0, (ii) when 
a specific relation (E = K for Uq and \M\ = \[K for Vq) exists between two constants 
of motion, and (iii) when the constant of motion £ 3 vanishes. In addition, all bounded 
trajectories are quasi-periodic rather than periodic in general and become truly periodic 
if a constraint condition on £ s (Eq. (2.34)), the same for Uq and Vq, is satisfied. In fact, 
there is an infinity of periodic trajectories in the neighborhood of a given quasi-periodic 
trajectory since the set of the rational numbers is dense in the set of real numbers. 

From a quantum mechanical point of view, the Uq system leads to a discrete spectrum 
while the spectrum for the Vq system comprises a discrete part, a zero point and a contin- 
uum. The Hamilton- Jacobi approach developed in sections 2 and 3, as a complement to the 
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earlier study of Kibler and Winternitz [20] , has allowed us to re-derive in a semi-classical 
way the discrete spectra for Uq and Vq. For both potentials, the accidental degeneracies 
corresponding to fixed values of \M\ and E are described by a unique dynamical invariance 
algebra, isomorphic to su(2), as shown by Quesne [2] for the Uq system and by Kibler and 
Winternitz [7] for the Vq system. This result is quite remarkable, especially in view of the 
fact that, in the limiting case (5 = 0, the dynamical invariance algebra is isomorphic to 
su(3) for Uq = o and to so(4) for Vq = q. 

As was mentioned by Kibler and Winternitz [7], the Vq system may for some energy 
levels exhibit higher degeneracies than those explained by the SU(2) dynamical invariance 
group. These are levels characterized by the triplets (m,n r ,ng) and (m 1 ,n' r ,n' e ), with 
m' 7^ ±m, satisfying 

4Q = - ( m + m') 2 ][I 2 - (m-m') 2 ] I = n' r + n' e - n r - n e ^ 0. (4.1) 

A similar result holds for the Uq system with / = 2n' p + n' z — 2n p — n z . 

Such degeneracies are "local" in that they are restricted to part of the energy spec- 
trum of the system (cf. Ref. [24]). The quantum mechanical operators X A related to a 
"local" symmetry would not commute with the Hamiltonian but would satisfy commuta- 
tion relations of the type 

[X A ,H] = \(H-E A ) (4.2) 

where E A are the energies of the levels for which the accidental degeneracy is observed. 
Since /, m, and m' are integers, Eq. (4.1) clearly imposes a constraint on the values of the 
coupling constant Q for which additional accidental degeneracy can occur. This type of 
constraint is reminiscent of Eqs. (2.36) and (3.31) that also impose conditions on Q. 

As pointed out recently by Moshinsky et al. [25], the problem of finding a group 
theoretical explanation of accidental degeneracy does not have an algorithmic solution. It 
is interesting to speculate about the relation between the periodicity condition (2.34) and 
the local degeneracies and we plan to return to this problem in the future. 
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Appendix : Torsion and Curvature for Hamiltonian Systems 



The purpose of this appendix is to derive closed formulas, in Cartesian coordinates, 
for the torsion and the curvature of the trajectories of a one-particle hamiltonian system. 
The formulas will involve solely the coordinates r(x, y, z) and the velocities r(x, y, z) of the 
particle. 

We start from the well-known formulas [26] 

(rAr).r |rAr 



|rAr| 2 | r* | 3 

for the torsion r and the curvature k of a curve (in R 3 ) in an arbitrary parametrization 
r = r(t). In the case where this curve is the trajectory of a particle (of mass //), we can 
use the Hamilton equations 

. OH . OH 

in order to express r and r in (Al) in terms of the Hamilton function H = H(jp, g, t). We 
shall employ the variable p and q in the Cartesian form (p x = \ix,p y = \iy,p z = ^z) and 
(x,y,z), respectively. Then, we have 

ldH ldH ldH 

x= — -K- y = — -7— z = — — (A3) 
H ox n oy n oz 

and 

1 fd 2 H d 2 H d 2 H d 2 H dH d 2 H OH d 2 H OH d 2 H\ 
x = 9 x + — — y + z — — — - — — — - — — — - — h 



\i \ dx 2 dydx dzdx dp x dx dx dp y dx dy dp z dx dz dtdx J 

(A4) 

with similar expressions for y and "z. As a result, the quantities r and k can be developed 
in terms of (x, y, z) and of the derivatives of H which turn out to be functions of (x, y, z), 
(x, y, z) and t. This leads to formulas, not reported here, which are intricate but very easy 
to handle with a symbolic computer language like Reduce or Macsyma. 

We now restrict ourselves to the case of conservative systems for which 

H = — {p x 2 + p y 2 + p z 2 ) +V(x,y,z). (A5) 
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Then, the expressions for 



— S den = (i 2 + » 2 + z 2 ) 3 ' 2 (M) 

can be given a simple form. Indeed, we get 

\i l x NUM = y — z— ~k~tX + -^-^-y + ^~^~z 



(A7) 



dz dy J \ dx 2 dydx dzdx 
dV dV\ ( d 2 V d 2 V d 2 V 



4 dx dz J \dxdy dy 2 ~^ dzdy 

(x— - ■— \ ( x + • + — i 

V <9j/ ^ dx J \dxdz X dydz^ dz 2 Z 

2 ^ AT f.dV .dV\ 2 f.dV .dV\ 2 ( . dV .dV\ 2 
/U 2 xDEN = [y— -z—\ + (z—-x—) + (x— - y—\ 

\ dz dy J \ dx dz J \ dy dx J 

or in compact form 

1 f dV dV dV \ 1 o 

NUM= — ~ (r A VV) . r-V— i+ r.V— j+ r.V— )k DEN = r A W| 2 

V ^ #2/ dz J [i z 

(A8) 

so that r and are given by closed form expressions involving only (x,y,z) and (x,y, z). 
Equations (A6) and (A7) are applied to V = Uq and V = Vq in the main body of this 
paper. 
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